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Abstract 

Perfect Space-Time Codes (STC) are optimal codes in their original construction for Multiple 
Input Multiple Output (MIMO) systems. Based on Cyclic Division Algebras (CDA), they are full-rate, 
full-diversity codes, have Non- Vanishing Determinants (NVD) and hence achieve Diversity-Multiplexing 
Tradeoff (DMT). In addition, these codes have led to optimal distributed space-time codes when applied 
in cooperative networks under the assumption of perfect synchronization between relays. However, 
they loose their diversity when delays are introduced and thus are not delay-tolerant. In this paper, 
using the cyclic division algebras of perfect codes, we construct new codes that maintain the same 
properties as perfect codes in the synchronous case. Moreover, these codes preserve their full-diversity 
in asynchronous transmission. 

Index Terms 

Cooperative Communication, Distributed Space-Time Codes, Perfect Codes, Delay-Tolerance, Cyclic 
Division Algebra, Tensor product. 

I. INTRODUCTION AND PROBLEM STATEMENT 

During the past decade, MEMO techniques have experienced a great interest in wireless 
communication systems. Using multiple antennas at the transmitter and the receiver provides 
high data rates and exploits the spatial diversity in order to fight channel fadings and hence 
improve the link reliability. Lately, cooperative diversity has emerged as a new form of spatial 
diversity via cooperation of multiple users in the wireless system [1J. While preserving the same 
MIMO benefits, it counteracts the need of incorporating many antennas into a single terminal, 
especially in cellular systems and ad-hoc sensor networks, where it can be impractical for a 
mobile unit to carry multiple antennas due to its size, power and cost limitations. 
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In cooperative networks, users communicate cooperatively to transmit their information by 
using distributed antennas belonging to other independent terminals. This way, a virtual MIMO 
scheme is created, where a transmitter is also acting as a relay terminal, with or without some 
processing, assisting another transmitter to convey its messages to a destination. The cooperative 
schemes have been widely investigated by analyzing their performance through different cooper- 
ative protocols Q]|-|[3]|. These protocols fall essentially into two families: Amplify-and-Forward 
(AF) and Decode- and-Forward (DF). In order to achieve the cooperative diversity, space-time 
coding techniques of MIMO systems have also been applied yielding many designs of distributed 
space-time codes under the assumption of synchronized relay terminals [|2]|-[|4]|. 

However, this a priori condition on synchronization can be quite costly in terms of signaling 
and even hard to handle in relay networks (5]|, O. Unlike conventional MIMO transmitter, 
equipped with one antenna array using one local oscillator, distributed antennas are dispersed 
on different terminals, each one with its local oscillator. Thus, they are not sharing the same 
timing reference, resulting in an asynchronous cooperative transmission. 

On the other hand, in a synchronous transmission, the distributed STCs are constructed basi- 
cally according to the rank and determinant criteria [7] and hence aim at achieving full diversity. 
Note that the rows of the codeword matrix represent the different relay terminals (antennas). 
So, when asynchronicity is evoked, delays are introduced between transmitted symbols from 
different distributed antennas shifting the matrix rows. This matrix misalignment can cause rank 
deficiency of the space-time code, and thus performance degradation. 

Therefore, the codes previously designed are no more effective unless they tolerate asyn- 
chronicity. Furthermore, an efficient code design should satisfy the full-diversity order for any 
delay profile. This intends to guarantee full-rank codewords distance matrix i.e., its rank equal to 
the number of involved relays, hence leading to the so-called delay-tolerant distributed space-time 
codes 10. 

II. Delay-Tolerant Distributed Space-Time Codes 

The first designs of such codes were presented by Li and Xia O as full-diversity binary 
Space-Time Trellis Codes (STTC) based on the Hammons-El Gamal stacking construction, its 
generalization to Lu-Kumar multilevel space-time codes, and the extension of the latter codes 
for more diverse AM-PSK constellations (HI, BUl. Systematic construction including the shortest 
STTC with minimum constraint length was also proposed in [[Toll , as well as some delay-tolerant 
short binary Space-Time Block Codes (STBC) ifTTTl . Recently, Damen and Hammons extended 
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the Threaded Algebraic Space-Time (TAST) codes to asynchronous transmission lfj~2ll . The delay- 
tolerant TAST codes are based on three different thread structures where the threads are separated 
by using different algebraic or transcendental numbers that guarantee a non-zero determinant of 
the codewords distance matrix. An extension of this TAST framework to minimum delay length 
codes was considered in lfT3ll . 

Meanwhile, perfect space-time block codes that are optimal codes originally constructed for 
MIMO systems lfT4ll - [fT7l . were also investigated for wireless relay networks. In lfi~8ll . lfi~9ll . the 
authors provided optimal coding schemes in the sense of DMT tradeoff EOl . based cyclic division 
algebras for any number of users and for different cooperative strategies. Nevertheless, all these 
schemes assumed perfect synchronization between users. Then, it was in [21J that Petros and 
Kumar discussed the delay-tolerant version of the optimal perfect code variants for asynchronous 
transmission. They stated that delay-tolerant diagonally-restricted CDA codes and delay-tolerant 
full-rate CDA codes can be obtained from previous designs by multiplying the codeword matrix 
by a random unitary matrix. This matrix can be taken specifically from an infinite set of unitary 
matrices that do not have elements in the code field. 

In this paper, we construct delay-tolerant distributed codes based on the perfect codes algebras 
from a different point of view. The new construction is obtained from the tensor product 
of two number fields, one of them being the field used for the perfect code. The codes are 
designed in such a way to maintain the same properties of their corresponding perfect codes 
in the synchronous transmission, namely full-rate, full-diversity and non-vanishing minimum 
determinant. In addition, unlike the perfect codes, the new codes preserve the full diversity in 
the asynchronous transmission. 

III. Background 

Before addressing the STC construction, we dedicate this section to briefly review the remark- 
able properties of the perfect codes as analyzed in [fT4ll - lfT7l . Then, following the framework of 
1161 , we present the cooperative communication model of interest. 

A. Perfect Space-Time Block Codes 

The concept of Perfect Code was originally proposed in lfl"4ll . [fT51l for N t = 2, 3, 4, 6 transmit 
antennas to describe a square N t x N t linear dispersion STC C. The perfect codes are constructed 
from cyclic division algebras ^4(K/F, a, 7) of degree n = N t defined by 
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- K and F are number fields and Ok, Of the corresponding ring of integers. F is called the 
base field and taken as F = Q(i) or F = Q(j) since the ST code transmits g-QAM or q- 
HEX information symbols for Nt = 2, 4 or N t = 3, 6, respectively. Thus, the constellations 
can be seen as finite subsets of the ring of Gaussian integers Of = Z[i] or Eisenstein 
integers Of = Z[j] [i = y/^T,j = e 27 ™/ 3 ), respectively. 

- K/F is a cyclic Galois extension of F of degree [K : F] = n with K = Q(i,9) or 
K = Q(j, 6) a field extension appropriately chosen in order to get an existing lattice and 
a division algebra, and 9 an algebraic number. 

- a is the generator of the Galois group Gal(K/F), Gal(K/F) = (a) = {a k }t =v For an 
element x G K, the conjugates of x are a k (x). So, the norm Nk/f and the trace Ttk/f are 
defined respectively as 

n n 

N K/¥ (x) =Y[<r k (x), Tr K/F (x) = 5> fe (x). (1) 

k=l k=l 

- 7 G F* = F/{0} the set of non-zero elements of F. It is a non-norm element suitable for 
the cyclic extension K/F lfi~5ll . 

The cyclic division algebra is then expressed as a right K-space 

A 1 = K © eK © e 2 K © . . . © e^K (2) 
with e e A, e n = 7 e K, 7 ^ and Ae = ecr(A) for all AgK. (3) 
The Perfect Codes V satisfy the criteria: 

• Full-rate: The code transmits iV 2 symbols drawn from QAM or HEX constellation and 
thus has a rate of R = N t symbols per channel use (spcu). 

• Full-diversity: According to the rank criterion Q, the determinant of the codeword 
distance matrix A = (Xj — Xj)(Xj — Xj)t for any two distinct codewords is non-zero. 
By code linearity, it can be reduced to 

det(A) ^ => det(XX f ) = | det(X)| 2 ^ 0, X ^ 0, X e C (4) 

• Non-vanishing minimum determinant: The minimum determinant of any codeword 
distance matrix, prior to SNR normalization, is lower bounded by a constant tp that is 
independent of the constellation size 

5(C) = min I det(X)| 2 > ^ > (5) 
o^xec 

• Cubic shaping: The QAM or HEX constellations are normalized according to the power at 
the transmitter so that the real vectorized codeword vectors are isomorphic to cubic lattices 
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Z * or A * . In other words, the rotation matrix M encoding the information symbols 
into each layer is required to be unitary to guarantee the energy efficiency of the codes. The 
shaping constraint leads thus to two other properties. The first one is the Uniform average 
transmitted energy per antenna. The second one is the Information losslessness as the 
unitary linear dispersion matrix M allows to preserve the mutual information of the MIMO 
channel. 

Thanks to prominent results on diversity-multiplexing tradeoff Il20ll . the perfect codes also 
verify two other equivalent properties: 

• DMT optimality: In [16], Elia et al. proved that the full-rate STCs from cyclic division 
algebra having NVD property achieve the optimal DMT in Rayleigh fading channel. 

• Approximate universality: Being CDA-based codes with NVD property, the perfect codes 
are approximately universal and achieve DMT for arbitrary channel fading distribution. 

Satisfying all these criteria, the perfect codes showed to improve the performance in terms of 
error probability upon the best known codes. 

B. Cooperative System Model 

In the sequel, we consider a cooperative system with a source S communicating to a destination 
D via M relays Ri in two phases as in Figure [Q and without direct links between the source and 
the destination. In the first phase, the source broadcasts its message to the potential relays. In the 
second phase, the relays use the DF protocol to detect the source message then if successfully 
detected transmit it to the destination. We assume that all the M relays are able to achieve 
error free decoding which could be possible by selecting the source-relays links, and consider 
only the links that are not in outage. Note that it could also be possible that not all the relays 
may successfully decode the original message, so the number of transmitting relays is usually 
assumed as a random variable. Since the relays transmission overlap in time and frequency, they 
can cooperatively implement a distributed space-time code. 

Considering only the second phase of transmission, the system is equivalent to a MIMO 
scheme where the distributed M x T perfect space-time code V is used by the relays, with 
M transmit antennas one by relay, and N r receive antennas at the destination. Every time slot 
t,t = 1 . . . T, the relays send the M x 1 t th column vector X t of the codeword X and the 
destination receives 

Y t = HA + W t , Y t , W t G C^ xl (6) 
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where Wt is the additive white Gaussian noise with i.i.d complex Gaussian variables with zero- 
mean and variance N , ~ A/" c (0, iV ), iV" = 2a 2 , a 2 being the noise variance per real dimension. 
H( represents the iV r x M complex channel matrix modeled as i.i.d Gaussian random variables 
with zero mean and unit variance ~ A/" c (0, 1). The channel is assumed quasi-static with constant 
fadings during a transmitted codeword and independent fadings between subsequent codewords. 
Dealing with square STCs (M = T), the codeword matrix X* contains M 2 information symbols 
si, . . . ,sm 2 carved from two-dimensional QAM or HEX finite constellations denoted by S. 

C. Asynchronous Cooperative Diversity 

The above expression © is valid only when relays are synchronized. In the presence of 
asynchronicity, the codeword transmission is spanned on more than T symbol intervals due 
to delays. Although the symbol synchronization is not required, we assume that the relays are 
synchronized at the frame-codeword level, which can be provided by means of network feedback 
signaling from the destination. Therefore, the start and the end of each codeword are aligned 
for different relays by transmitting zero symbols, and hence there is no interference between 
codewords transmission. We further assume that the timing errors between different relays are 
integer multiples of the symbol duration and the fractional timing errors are absorbed in the 
channel dispersion. In the codeword matrix, these delays are also filled with zeros; they are 
known at the receiver but not at the transmitting relays O. 

Denoting a delay profile by D = (Oi, D 2 , . . . ,0m), a delay 0; corresponds to the relative delay 
of the received signal from the i th relay as referenced to the earliest received relay signal. 
Let o max denotes the maximum of the relative delays, then from the receiver perspective, the 
M x (T + max ) codeword matrix was sent instead of the M x T space-time code. 

D. Motivation of the Code Construction 

The diversity order of any space-time code is defined by the minimum rank of the distance 
codeword matrix over all pairs of distinct codewords [7|. The distributed M x T perfect codes V 
are full-rate full-diversity for the synchronous transmission between the relays and the destination. 
Note that in general, a transmission between source, half-duplex relays and destination will result 
in rate loss. When asynchronicity is introduced, the code is no more full-rate since it is spanned 
on (T + max ) time instants. Moreover, certain delay profiles D can result in linearly dependent 
rows, thus the code will loose its full-diversity property. Let us illustrate this by the following 
example. 
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(7) 



Example of Golden Code: We consider the distributed 2x2 Golden code transmitting 4 
information QAM symbols si, s 2 , s 3 , s 4 from two synchronized relays with the codeword matrix. 

+ s 2 9) a(s 3 + s A 9) 
ia(s 3 + sj)) a(si + s 2 6) 
The Golden code is designed on a cyclic field extension of degree 2 over the base field Q(i). 
Using the generator matrix of the corresponding complex 2-dimensional lattice, the codeword 
elements are lattice points obtained by linear combination of pairs of symbols. 

Now, let the first relay be delayed by one symbol period with respect to the second d = (1,0), 
such that the new asynchronous codeword matrix be 



Xa "7! 



(8) 



a(si + s 2 6) a(s 3 + sj)) 

ia(s 3 + sj)) a(si + s 2 6) 
Suppose we have two distinct codewords Xi and X 2 with s^i = — si >2 = — si and the other 
symbols equal i.e., s i:1 = s ij2 ,i = 2 . . .4. The difference between matrix codewords is defined 
in both synchronous and asynchronous cases as 



A(s), = 



2asi 
2as 1 



A(s) a = 



2asi 
2as 1 



(9) 



It can be seen that A(s) s is a full-rank matrix whereas A(s) a has rank one, so the Golden code 
is not a delay-tolerant code. 

In fact, it can be seen from the asynchronous codeword matrix X a that some symbols are 
aligned at the same instant due to delays loosing thus diversity. In order to resolve this problem 
of rank deficiency, our solution consists in transmitting from each antenna (relay) at each 
transmission time a different combination of all the 4 information symbols. This way, in the 
presence of delays, we ensure that any combined symbol sent from the 2 relays arrives at the 
destination in at least 2 different instants, hence guaranteeing the full-diversity order of the 
space-time code. 

A new 2x2 STC will have then the shifted codeword matrix 

fi(s 1 ,s 2 ,s 3 ,s A ) / 2 (si, s 2 , s 3 , s 4 ) 

73(51,52,53,54) h(s 1 ,s 2 ,s 3 ,s 4 ) 
Now, to get these 4 linear combinations of the 4 symbols, we need a higher dimensional lattice 
(n = 4) compared to the 2-dimensional lattice used for the Golden code. So, we propose to obtain 
the corresponding 4x4 lattice generator matrix by the tensor product of two field extensions of 
Q(i), one of them being the field extension of the Golden code. 



(10) 
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Following this idea, we aim at constructing, in general, new M x M codes that are based on 
CDA of the M x M perfect codes such that they maintain the same optimal properties as perfect 
codes in the synchronous case. But also, these codes preserve their full-diversity in asynchronous 
transmission and thus are delay-tolerant for arbitrary delay profile. 

IV. Construction of Delay-Tolerant Distributed Codes 
based Perfect Codes Algebras 

A. General Construction 

The approach consists in constructing a division algebra isomorphic to the tensor product 
(also called Kronecker product or cross-product) of two number fields of lower degrees. Other 
constructions based on the crossed-product algebras have been investigated in G2ll . ll23Tl either 
for prime or coprime degrees of the composite algebras. In these constructions, the space-time 
code was built on the cyclic product algebra. However, in the present construction, the higher 
degree algebra is only used to derive appropriately the space-time code. 

Since we intend to construct a full-rate M x M space-time code that is based on the CDA of 
the full-rate M x M perfect code, then the first algebra to be considered is the cyclic division 
algebra of the perfect code At(Ki/F, cr 1 , 7^ of degree n\ = M over the base field F. For sake 
of simplicity, we analyze in the sequel the case of Gaussian Field F = Q(z) to explain the 
construction. Indeed, we consider the cyclic field extension Ki = Q(i, 61) of degree n\ = M 
over F, 9\ being an algebraic number. The principal ideal X Kl C Kl is generated by an 
element a and its integral basis is Bi = (t>i,t>2, • • • ,%) (or if unitary, it is given by Bi = 
(a, aOi, . . . , a^ -1 )). The basis of the complex algebraic lattice A(Zk 1 ) is obtained by applying 
the canonical embedding to B^ Consequently, the generator matrix corresponds to the rotation 
matrix in Z[i] M 

vm 
Ci(vm) 

(11) 

v?-\v M ) _ 

where y/pi is a normalization factor used to guarantee the matrix unitarity. 

Now, we consider another Galois extension K 2 over F of the same degree n 2 = M such that its 
discriminant is coprime to the one of Ki i.e., (g?Kd ^k 2 ) = 1- L et K 2 = Q(9 2 ) with 9 2 an algebraic 
number. The Galois group is generated by a 2 as Gal(K 2 /F) = (a 2 ). The principal ideal of the 



Mi 



1 

Pi 



Vl v 2 
t7i(Vl) <Ti(v 3 ) 

1 [Vl) 1 (V2) 
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;ia, 



M-X\ 



algebra is such that Xk 2 = Ok 2 an d thus its integral basis is given by B 2 : 
The canonical embedding of B 2 gives another complex rotated lattice of Z[i] M that is generated 
by the unitary matrix M 2 with Jp~2 the normalization factor, 



M 2 



02 



1 <- 1 (9 2 ) 



On 



(12) 



The tensor product of both field extensions allows to build a rotated lattice in higher dimension 
corresponding to the complex M 2 x M 2 unitary matrix M based on the previous M x M 
constructions. According to Il24ll . 

Proposition 1; : Let IK be the compositum of the above Galois extensions, K = K]K 2 = 
Q(i, 9i, 62) of order n = nin 2 = M 2 over F as presented in Figure [2J 

Since Ki and K 2 have coprime discriminants, the corresponding lattice generator matrix can 
be obtained as the tensor product of the previous unitary generator matrices. 

1 



M = M 2 ® Mi 



VP1P2 



a 



M-l. 
1 1 



«1 



a 



.M-li 
1 



[Vi 



VM 



a 



M-X, 
1 1 



v M ) 



vm 



(Jo 



M-X, 



{vm) 



-M-X/., \nl 
1 (^l)^2 



„ -M-X/nM-X\ 

vio 2 [p 2 ) 



a 



1 



[vi)a 2 



M-X/qM-X^ 



0, qM-x 

Vm&2 



{v M )0{ 



., -M-X/nM-X^ 
VM&2 \?2 



a 



.M-X, 
1 



[v M )a 2 



M-X/nM-Xs 



(13) 



Consequently, 

Proposition 2: : Let rrij = [K : Kj] = n/rij, j = 1, 2 the order of the extensions, then the 
discriminant of K is g?k = d^d^\ The minimum product distance of the lattice is derived from 
the discriminant of K as 



d 



1 



1 



p,mm 



11112 

l K 2 



(14) 



[~T~ I 'jmi jr> 

Using the matrix M, the space-time coded components are given by the linear combination 



Ai 2 



x = Ms where s = [si, s 2 , . . . , sm 2 ] t is the information symbol vector carved from a g-QAM 
constellation (e Z[i] M2 ). Then, the space-time codeword matrix is defined by distributing the 
components with appropriate constant factors (pi, I = 1,...,M 2 . It can be represented as a 
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Hadamard product 

4>iX\ <f>M+l%M+l ■ ■ ■ 4>M(M-l)+lXM{M-l)+l 

- , 02^2 0M+2^M+2 4>M(M-l)+2 x M(M-l)+2 

X = [$] • [x] 



(15) 



4>M%M 4>2MX 2 M ■ ■ ■ 4>M^X M 2 

The key idea in the code construction is to determine the coefficients (pi that allow one to 
preserve the same properties of the corresponding perfect codes in synchronous transmission 
(Section UTTAl) . 

• On one side, it can be seen that the new code transmits M 2 information symbols and thus 
is full-rate with R = M spcu for a relays-destination transmission phase. 

• On the other side, we need to find the (pi factors that satisfy the rank criterion © in order 
to have full-diversity codes. 

• Moreover, the perfect codes have non-vanishing minimum determinants. Then, we are 
interested in deriving MxMST codes that have not only non-zero determinants, but also 
these determinants do not vanish when constellation size increases. 

• In order to guarantee uniform energy distribution in the codeword, we ask that (pi verify 
|0; | = 1. Choosing further the coefficients (pi E Of = Z[z] yields better determinants as 
obtained for the non-norm elements 7 of the perfect codes lfl5ll . This restricts the values 
of (pi to (pi = ±1, ±i. 

• It can also be noticed that the new code satisfies the cubic shaping property since the 
generator matrix M of the M 2 -dimensional lattice is unitary, and hence the code is 
information lossless. 

In addition, when asynchronicity between relays is involved, the rank criterion should be also 
verified for the shifted matrix and another criterion will be analyzed that is the non-zero product 
distance of the codeword matrix in order to prove that the new codes are delay-tolerant, and 
thus keep their full-diversity in asynchronous transmission. 

V. New Delay-Tolerant Codes from 2, 3, 4-dimensional Perfect Codes 

Based on the previous approach, we consider the perfect codes proposed in lfl4l . lfT51 for 
dimensions M = 2, 3, 4 to construct the new delay-tolerant codes. Then, in the next section, we 
apply this construction for the perfect codes presented for any number of antennas in [?]. 
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A. 2 x 2 Code based on Golden Code 

The Golden Code was constructed in [14J using the cyclic division algebra Ai(Ki/Q(z), oi, i) 
of degree 2 over F = Q(z). Ki = Q(i, a/5) is a Galois extension of degree 2. It is a 2-dimensional 
vector space of Q(i) with basis B = (l,#i), Q\ = 1+ 2 being the Golden number. Its Galois 
group Gal(Ki/F) is generated by o~i : v5 >-)■ — \/5. In order to get a rotated lattice A(Xk 1 ) 
of Z[i] 2 , the principal ideal Z Kl = aCid generated by a = 1 + i — iQ\ was found. Its basis is 
E>! = (a,a9i) and its unitary generator matrix is given by 

a a6i 



l-\/5 



and a = 1 



i6i the respective conjugates of 6i and a. 



(16) 



with 6i — 

Let K 2 — Q(^) the cyclotomic extension of degree 2 over F with 9 2 = Cs — ei7r ^ 4 the 
primitive 8 <h root of unity. Its discriminant d^ 2 = 4 and it is coprime to the one of Ki since 
c?Kj = 5. The Galois group Gal(K 2 /F) is generated by a 2 : Cs ^ — Cs an d the integral basis of 
K 2 is B 2 = (1, Cs)- The corresponding unitary generator matrix is 



1 Cs 
1 -Cs 



(17) 



Therefore, K = KiK 2 = Q(z,6 , i,Cs) is the compositum of Galois extensions of degree 2 
each, with coprime discriminants. The 4x4 unitary matrix is obtained by the tensor product of 
previous matrices as 

a aOi «Cs 



M 4 



1 



a^iCs 
a^iCs 



a aOi — «Cs — odiCs 
a aOi — aCs — Q^iCs 



(18) 



and the codeword matrix is defined by 

r(s) = r ^, (19) 

2 X 2 4 X 4 

where Xj are the components of the vector M 4 s with s 1; s 2 , s 3 , s 4 are g-QAM symbols. We pro- 
pose now to determine the coefficients (j>i, I — 1, . . . , 4 that satisfy the non-vanishing determinant 
criterion. 
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1 ) Non-vanishing minimum determinant: The determinant of this codeword matrix equals 



8(s) = (j)i(() A XiX A - 0203^2^3- 



By developing x\x± and x 2 x 3 , we obtain 

1 



1 + i VEa ia <- 



v/2 



2^3 



1 / 1 + i 1 + i 

— N Kl/¥ (a) \ G(s) -^=r\/5s 1 S 4 + —^=-^5233 



(20) 

(21) 
(22) 



with 



N^_ 1 /^(a) = aa = 2 + i and G(s) = s\ — s\ — is\ + is\ + S1S2 — is^s^. 



(23) 



It is interesting to note that the Golden codeword given by matrix CO) has a determinant of 



S'(s) = -N Kl/¥ (a)G(s). 



(24) 



Therefore, by choosing <pi = 3 = 4 = 1 and </> 2 = —1, the determinant of the new code 
is equal to the Golden code determinant, and does not vanish when increasing the size of the 
QAM constellation carved from Z[z]. Hence, the new code achieves the diversity-multiplexing 
tradeoff |fl6l . 

It can also be noticed that the coefficients <pi can be changed equivalently to the coefficients of 
the Fourier matrix F n = {w^ k ) where w = e 2l7r// " is the primitive n th root of unity. For dimension 
2, we have 

' 1 1 

1 -1 

Furthermore, we have find fixed unitary matrices U and V such that V = UGV for all values 

of si, s 2 , s 3 , s 4 with 



(25) 



U 



Cs o 






-1 



1 -1 



(26) 



2) Delay-tolerance: In the distributed setup, each row of the code matrix is transmitted by a 
different relay (Section ttll-BI) . In practical scenarios, the two relays do not share a common timing 
reference, and therefore, the arrival of packets is not synchronous. As we assume synchronization 
at the symbol level, the distributed code can still achieve full diversity if the differences between 
matrix codewords are full rank even when the different rows are arbitrarily shifted. In what 
follows, we prove that the new code T satisfy this condition. 
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(27) 



Consider the shifted codeword matrix of T 

X\ x 3 

—%2 X4 

we need to guarantee that it is full rank when s ^ i.e., for any s x 7^ s 2 from the constellation 
(rank = min(M, T + 1) = 2). This restricts to show that the 2x2 submatrix 

x 3 
-x 2 

is full rank i.e., its determinant x 2 x 3 7^ when s 7^ 0. 

More generally, having delay profiles D = (1,0) or (0,1), the problem turns to prove that 
the product distance in the rotated constellation associated with the matrix M 4 of T is non-zero 
over Z[i], so that any component product is non-zero. This product distance is evaluated as 



d r . 



n 

j=0 



\X 1 X 2 X 3 X4\ 



10 



-aa 



-aa 



1 + I r- 1 + { r- \ 

-V0S1S4 H —Vos 2 s 3 ) 



(G(s 
1 

" 20 



V2 



s/2 



x 



V2 



V5( 



S1S4 — S2S 3j 



10 --v-v~/ ^ V2 

with G(s) = sf — si — is\ + is\ + sis 2 — IS3S4 G Z[i] for s G Z[i] 4 . 

As a direct consequence from the tensor product construction, Equation (TT41) gives 



(28) 



p,min 



1 

20 



Thus, the minimum product distance is non-zero. It can also be verified in d p by setting si = 
1) s 2 = S3 = S4 = 0. So, dp is non-zero unless Si = s 2 = S3 = S4 = 0, and consequently the 
submatrix is full rank since x 2 x 3 7^ unless s = 0. 

Therefore, the new code unlike the Golden code keeps its full-diversity in the case of asyn- 
chronous relays. However, we cannot guarantee the non- vanishing determinant property in the 
asynchronous case because the expression of x 2 x 3 can be interpreted as a Diophantine approx- 
imation of by rational numbers which can be made tighter by using larger constellation 
size. 



B. 3 x 3 Code based on 3 x 3 Perfect Code 

In order to construct the delay-tolerant 3x3 code, we consider the base field F = Q(j) 
and we use g-HEX symbols. Let 6\ = (7 + Cf 1 = 2cos(y : ), with (7 the 7 th root of unity. 
The 3x3 perfect code was constructed using the cyclic division algebra ^4 1 (K 1 /F, j) of 
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i,0i) and <j\ the generator of the cyclic 



k=l 



order 3 lfi~5l . where the relative extension K% = 

extension Ki/F with a 1 : £ 7 + Cf 1 ^ C7 + Cf 2 ■ The integral basis is given by Bi = {t> fc } 
{(1 +j) + 9 1 , (-1 -2j) +j6l (-1 -2j) + + (1 +j)6l} and the complex lattice A{X Kl ] 

is a rotated version of Z[j] 3 . It is generated by 



Mi 



1 

7f 



Ul v 2 v 3 
<?i{vi) (T\{v 2 ) cri(f 3 ) • ( 29 ) 
cr^Oi) of(> 2 ) cr?(w 3 ) 

The relative discriminant of Ki is d^ x = 49. Another extension of F of degree 3 that has 
coprime discriminant with Kq is the cyclotomic extension K 2 = Q(Cq) with C9 = e 2in ^ 9 the 
primitive 9 th root of unity and d K2 = 27. Its Galois group Gal(K 2 /F) is generated by cr 2 : C9 h> 
j'C9- The integral basis of K 2 is B 2 = (1, ( 9 , Cf ) and the lattice generator matrix is 



c 9 c, 



1 a, fCi 



(30) 



1 j%9 jQ 

The compositum of both extensions K = KiK 2 = Q(j, 2 cos(y-), C9) is of order 9 over 
Then, the corresponding 9-dimensional complex lattice is generated by the 9x9 unitary matrix 



Mc 



1 



c 9 c 



1 K* j 2 Q 



V-2 



v 3 



aAVx (Ti(v 2 crAv 3 



(31) 



1 f(9 jQ 

and the 3x3 space-time code is defined by the matrix 

01^1 04^4 07^7 

Tfs) = (j) 2 x 2 5 x 5 8 x 8 , (32) 

03^3 06^6 09^9 

where x, are the components of vector M 9 s, s being the information symbol vector carved from 
g-HEX 9 constellation. 

1 ) Non-vanishing minimum determinant: By proceeding as previously, we need to determine 
the coefficients 0/, / = 1, . . . , 9 that guarantee the non-vanishing minimum determinant. In order 
to get |0; I = 1 so that a uniform average energy is transmitted per antenna, and to obtain better 
values of the determinant, we limit the choice of 0; to 0/ = ±1, ±j. 

By developing the code determinant using symbolic computation under Mathematica, we find 
that it has the same expression as the 3x3 perfect code determinant by choosing 0/ as the 
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Fourier matrix coefficients in 



1 1 1 

<[, = 1 j j 2 

j 

Therefore, the 3 x 3 infinite code T(s) has non-vanishing minimum determinant equal to 



1 J 2 



(33) 



due, 49 



(34) 



2) Delay-tolerance: On the other hand, to prove the delay-tolerance of this code, we should 
guarantee that the corresponding shifted codeword matrices are full rank. Therefore, it suffices 
to verify that for each asynchronous matrix there exists a square 3x3 matrix that is full rank 
i.e., its determinant is non-zero. In fact, if we enumerate all the delay profiles, it can be noticed 
that the problem of guaranteeing full-rank shifted matrices turns to guarantee that 

- All component products C d p are non-zero. This condition is always verified since the 
product distance d p = lf i=1 \x { \ ^ over Z[j] as d p , min = V49 l 2?3 ■ 

- All 2x2 minors of T(s) are non-zero that is equivalent to verify that the 9 entries of the 
cofactor matrix of V are non-zero. 

In order to prove the second condition, we find two unitary matrices U and V such that the 
codeword matrix V can be written as V = UZV for all s, with Z is the perfect code matrix and 
U and V are defined by 



U = 







o fc? 









9 

o j% 



V 



1 

7! 



i i i 

Cg j'Cq J 2 Cq 

^9 J S9 J^9 



(35) 



Let define the cofactor matrix of the perfect code by Z. Since Z is a finite subset of the cyclic 
division algebra Ai, Z is also a subset of Ai taken from the lattice A = (9 Kl © e(9 Kl © e 2 OK 1 - 
with e 3 = j and Kl is the ring of integers of K x . Hence, the cofactor matrix can be represented 
as a 3 x 3 codeword matrix. For simplicity, we denote by z = cri(z) and z = af(z), the conjugates 
of an entry of the codeword matrix. The cofactor codeword matrix is then defined by 



Z = 



z l z 2 z 3 

jz 3 Zi z 2 
jh 3% zi 



(36) 
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where each diagonal Z, = Mi[sj, Sj+i, Si +2 \ T , i = 1, . . . , 3. 

Since T = UZV, we denote T its cofactor matrix. It is given by T = V^ZU^ and satisfies 



IT = UZVV+ZU 1 = det(Z)I, 



(37) 



with 



1 

o c 9 o 
o o cl 

Developing the cofactor matrix T, we get 



Vt = ^ 



1 






1 


Ki 




1 


a 





(38) 



r = v f zu f 



Z\ + Qz 3 + (9Z2 C9Z2 + Zi + Cg^3 CqZ 3 + ( 9 z 2 + z x 



Zi + j %9 Z 3 + jCg^2 Cg^2 + r^i + jCg^s Cg^s + 3%%^ + jz x 
zi + jQzz + j 2 (9h C9Z2 + jzi + fQh Qz z + j( 9 z 2 + J 2 Ii 



•(39) 



Note that the Galois group G = Gal(K/F) has two generators o\ and o 2 , it is given by 

G = Gal(Ki/F) x Gal(K 2 /F) = (o x , o 2 ) = {1, a u a 2 , of, a 2 , o x o 2 , o x o\, o\o 2 , o\o 2 2 }. (40) 
From the expression of T (1391) , we define 

X x = zi + C9Z3 + (9^2, (rio 2 (Xi) = zt+ fQh + 3(9*2, 0I0KX1) = z x + j 2 C 9 ^2 + j(gZ 3 , 
X 2 = z x + Cg^2 + Cg%, o 1 o 2 (X 2 )=z 1 + j( 9 z 2 + J 2 C|^3, (y\o\{X 2 ) = z x + j 2 C 9 I 2 + jCgZ 3 , 

X 3 = \ + C 9 ^2 + Qz 3 , (T 1 0r 2 (X 3 )=2l + j'Cg^2 + J 2 C 9 ^3, ^2(^3) = *1 + /Cu^fe + jCf^l) 

the elements and their conjugates by the embeddings [o\0 2 ) k , k = 0, . . . , 2 with a^cr^Xj) = 
o\o\{Xi) = oio 2 (Xi). We also have X 2 = oi(X x ) and X 3 = o\(Xi) the conjugates of X\ 
by the embeddings o\. Then, the cofactor matrix can be rewritten as 



X\ X2 X 3 

cricr 2 (X 3 ) j 2 oxo 2 (Xx) 30\o 2 {X 2 ) 
o\o\{X 2 ) joloi(X 3 ) fofoKXx) 



Xx ox(Xx) ofiXx) 
o 2 {Xx) j 2 oxo 2 (Xx) jo\o 2 {Xx) 
oKXx) joxoi(Xx) faloKXx] 



(42) 



Finally, computing the product distance of this matrix, we get the product of all the 9 entries 
Xi G K that is the product of all the 9 conjugates of X\ G K, and thus 



dp = 11 1^*1 = I^M^OI 



(43) 



,_, v^k V49 3 27 3 

As a result, the elements of T are all non-zero unless s = which concludes our proof on the 
full-diversity of the 3 x 3 code T, hence its delay-tolerance for any arbitrary delay profile. 
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C. 4 x 4 Code based on Ax A Perfect Code 

Similarly to the 2x2 case, the 4x4 code is derived over F 



(i) based on the 4 x 4 
,,,, ^ , ^ - Q(i,0O 

of degree [Ki : F] = 4 and its relative discriminant is d Kl = 1125. The cyclic Galois group 



perfect code algebra. Let 9i = C15 + C15 1 — 2cos(ff), the relative extension is K4 



k=i 



Ki/F is generated by o\ : C15 + C15 1 ^ C15 + Ci5 2 - The integral basis is Bi = {v k } 
{(1 - 3i) + i0f , (1 - 3i)0i + -i + (-3 + 4i)0i + (1 - i)0?, (-1 + i) - 30i + 0? + 0?} and 
the complex rotated lattice of Z[i] 4 is generated by the unitary matrix 



Mi 



(44) 



Vl v 2 v 3 v 4 

1 cti(di) <7i(i> 2 ) <ti(v 3 ) ai(v 4 ) 

<T?(ui) <T?(U 2 ) of^i) <Ti(«4) 
<T?(ui) <T?(U 2 ) <7?(u 3 ) <T?(U 4 ) 

The second relative extension K 2 is chosen such that its degree is 4 over F and has coprime 
discriminant with K4. Let K 2 = Q(Ci6) this cyclotomic extension with d^ 2 = 256 and Ci6 = e l7r / 8 
the primitive lQ th root of unity. The cyclic Galois group is generated by a 2 : Ci6 ^ «Ci6- The 
integral basis of K 2 is B 2 = (1, (f 6 ) and the lattice generator matrix in Z[i] 4 is given by 



M, 



Cie 

<16 
-Cl6 
-<16 



t 2 
S16 

S16 

S16 



S16 



16 



-c 3 



16 



-Ci 2 6 K 



16 



(45) 



Then, the tensor product of both cyclic extensions defines the compositum field K = KiK 2 = 
Q(i, 2cos(||), Ci 6 ) of order 16 over Q(i). Accordingly, the 16-dimensional complex lattice is 
generated by the 16 x 16 unitary matrix M 16 = M 2 <g)Mi. The 16 codeword elements are derived 
from the linear combination M i6 s of q-QAM information symbols. They are then distributed in 
the 4 x 4 codeword matrix and assigned the coefficients 0^, I — 1, . . . , 16 as 

(f) 1 X 1 5 X 5 (f) 9 Xg 0i 3 Xi3 

(f) 2 X 2 6 X 6 010^10 014^14 

3 X 3 7 X 7 011X11 015^15 

4 X 4 08X 8 0l 2 X i2 0l 6 Xi6 



r(s) 



(46) 



1) Non-vanishing minimum determinant: The coefficients <pi are restricted to \<pi\ = 1 for 
uniform energy transmission and should satisfy the NVD criterion. Therefore, as in previous 
dimensions, computing the code determinant using symbolic computation under Mathematica, 
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we find that such coefficients corresponding to the Fourier matrix coefficients in Q(i) allow to 
get a 4 x 4 space-time code with the same determinant as the perfect code. We have 

1111 

1 i -1 -i 
1-11-1 
1 -i -1 i 

Therefore, the 4x4 infinite code T(s) has non- vanishing minimum determinant 

1 1 



(47) 



^min(C) 



d Kl H25' 

and the 4 x 4 codeword matrix is defined for xi = X by 

■( ! .)•<, -r !:; X <r 2 (X) 

.1-, /.('(I -.r ;L o -/.C14 °"iP0 icr 2 cri(X) 

x 4 -ixg —x 12 ixie 



(48) 



r( s ) 



af(X) -a 2 al(X) 



ai(X) 
a\a\{X) 



a\{X) 
-io\o\{X) 



3 2 

-(To a 



1(X) 



af(X) -i<r 2 oi(X) -cr 2 2 af(X) io%of(X) 



(49) 



2) Delay-tolerance: Now, let us examine the delay-tolerance aspect of this code. For this task, 
we start by enumerating all the types of delay profiles. Consider the integer numbers a, b, c, d 
with a b ^ c ^ d and < a, b, c, d < 3, we can define four types of profiles as: 

- Type 1 of form Z> = (a, 6, c, d) 

- Type 2 of form Z> = (a, a, 6, c) 

- Type 3 of form Z> = (a, a, b, b) 

- Type 4 of form Z> = (a, a, a, 6) 

Each of the asynchronous shifted codeword matrices corresponding to these profiles is full 
rank if and only if it includes a square 4x4 matrix that is full rank i.e., a 4 x 4 minor that is 
non-zero. This will be proved in the sequel for the different delay profile types. 

Types 1 and 4 

If we consider the delay profiles of types 1 and 4, for instance Z)i = (0, 1, 2, 3), ~d 2 = (0, 0, 0, 1) 
and Z) 3 = (3, 0, 0, 0), the 4 x 4 minors A^ 4X 4 relative to the 4 x 4 shifted matrices can have one 
of these expressions 

- The product of some components of the codeword matrix T: J j \4>iXi\ 

1<«<16 

- The product of one component and a 3 x 3 minor Ms xS : 4>iXiM.^'i 
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Proof 1 



16 



In the first case, we have by construction that all component products C d p = ] \ \x k \ are 

k=l 



non-zero since d Pjia i n = -k= 



i_ 



Vll25 4 256 4 ' 

In the second case, following the same analysis of the 3 x 3 space-time code, we find the 
unitary matrices U and V 



U = 



1 













l 


1 


1 


1 















Cl6 




-Cl6 


-<16 















c 2 

S16 


-C 2 

S16 


c 2 

S16 


-C 2 











-<16 . 




C 3 


-<1 3 6 


-C 3 
S16 


<1 3 6 



(50) 



such that the new 4x4 code T can be written as T = UZV, Z being the 4 x 4 perfect code. 
Then, we derive the cofactor matrix T and prove that it has non-zero entries as its product 
distance is non-zero. Thus, the 3 x 3 minors are full-rank yielding full-rank shifted matrices. 

Type 2 

For delay profiles of type 2, for instance Z)x = (0, 0, 1, 3), Z> 2 = (2, 2, 1, 0) and D 3 = (3, 3, 0, 2), 
we can find 4x4 minors A^ 4X 4 in the relative shifted codeword matrices that are equal to 



M 4x4 = {(t>jXi){<f>kXk)M 2 x2, 



(51) 



where M.2x2 has its components <pixi such that only one <pi = ±i. So, the 4 x 4 minors are 
non-zero if these 2x2 minors are non-zero for any s ^ 0. 

Proof 2 



Mi 



Let 

Xi x 5 
x 2 ix 6 

be such 2x2 minor and consider any 3x3 minor A^3 X ,3 that includes Mi, for example 

Xi x 5 x 9 

:>> :i - X 2 IX 6 — X w 
X 3 -X 7 Xu 

It can be expanded into 

x 5 x 9 
ix 6 -x 10 



(52) 



M 



M 3x3 = x n 


X\ 


X 5 




X\ 


Xg 








+ x 7 






+ x 3 




%2 


IXq 




x 2 


-x 10 





x u Mi + x 7 M 2 + x 3 M 3 



(53) 
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By developing the 2x2 minors, we have according to the 4x4 codeword matrix (l49l) 

Mi = iXa 2 ai{X) - a x {X)a 2 {X) 

M 2 = -Xal<Ji{X) - ai(X)a 2 2 (X) = -Y - a 2 2 (Y) 

M 3 = -a 2 {X)ala 1 {X)-ia 2 ai{X)al{X)=ia 2 {Mi) (54) 

then 



-M.3X3 = aia((X)Mi + a 2 af(X)M 2 + i<rt(X)<T 2 (Mi). 



(55) 



If Mi = 0, A^3x3 = o 2 o\[X)M 2 can be zero since M 2 is a trace and can be zero (if 
Y — 7Ci6)7 ^ ^i,Y G K). However, we have from Proof 1 that any 3x3 minor is non- 
zero over Z[z]. Thus, Mi cannot be zero over Z[z] unless s = 0. By a similar analysis, we can 
prove that any 2x2 minor of the same form of Mi is non-zero for s 0. 

Type 3 

For this type, we distinguish two cases of profiles: 
31- a = 2, 3, b = such as Di = (2, 2, 0, 0), Z> 2 = (3, 3, 0, 0) 
311- a = 1, b = such as D = (1, 1, 0, 0) 

In the first case, there exist 4x4 minors that are equal to the product of two 2x2 minors 
such that these M 2X 2 have there components 4>iXi with only one 4>i = ±z, hence are non-zero 
according to Proof 2. 

In the second case, the 4x4 minors are functions of 2 x 2 minors M 2x2 as following 



2 1 

-M 4X 4 = TY\M 



k=0 1=0 



2x2,fc+/+l) 



(56) 



and thus we have to prove that this sum is non-zero over Z[i}. For this task and without loss of 
generality, we consider the delay profile D = (1, 1, 0, 0). 
Proof 3 

Let the 4x4 minor relative to this delay profile be 








Xi 


x 5 


Xg 


M Ax4 = 





x 2 


ix Q 


— XiQ 




X3 


-x 7 


Xix 


-Xl5 






-ix s 


-Xl2 


ix 16 



Xi x 5 
x 2 ixQ 



X 3 ~Xl5 

X4 ixiQ 



Xi Xg 
%2 —XlO 



X3 Xu 
X4 —Xi 2 



X 5 Xg 
IX 6 -XiQ 



MiM 2 - M 3 M 4 + M 5 M 6 , 



x 3 -x 7 
x 4 — ix$ 

(57) 
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with according to the codeword matrix in Equation (l49l) 



then 



Mi = 


■ iXa 2 a 1 (X)-a 1 (X)a 2 {X) 


M 2 = -ialal{Mi) 




M 3 = 


-Xaia l (X)-a 1 (X)al(X) 


M 4 = al(M 3 ) 




M 5 = 


i<r 2 (Mx) 


M 6 = -al(Mi) 


(58) 


M ixi 


= -iM 1 alal(M 1 ) - M 3 a\{M 3 ) 


- ia 2 (Mi)al(Mi). 


(59) 



By denoting the first term in this expression V\ and the second term V 2 , then 

M AxA = V l +V 2 + a 2 (V 1 ). (60) 
Recalling M 3 , we can notice that it can be written as 

M 3 = -Y - a 2 2 (Y) = -Tt k/¥i (Y) with Y = Xa 2 2 a 1 (X) G K. 
Let Y G K be Y = A x + Cie^i + Cfed + Cw D i with A U C u £>i G K x . Then, 

A4 3 = -2Ax - 2Ci 6 Ci = ai + Cs&i with ai, &i G Ki. (61) 
For simplicity, we denote the conjugate a\(x) = x, so we have 

V 2 = -M 3 al(M 3 ) = -(^ + C 8 6i)(Si + (sh) = -aiaj - Cl&i&a - Csfai&i + (62) 

Let us now examine the nested sequences of fields included in the compositum field K in 
Figure [3] We have 

F = Q(z) c L x = Q(z, v/5)cK! = Q(z, 0i) CF 1 = Q(z, ^( 8 )cK = Q(z, 1; Cie) (63) 
F = Q(z) C L 2 = Q(z, (s) CK 2 = Q(z, Ci 6 ) C F 2 = Q(z, Cie, \/5) Cl = Q(z, X , Cie) (64) 

with the perfect algebra VA = (Ki/Li, of, 71) = Ki © W1K1, where u 2 = 71 = 2 and <xi : #i = 
C15 + Cis 1 ^ 9f - 2 = C 2 5 + Q 2 , \/5 1 — ^ — \/5. As we have oi(\/5) = -a/5, Li is the subfield 
fixed by (<x 2 ) the subgroup of order [L x : F] = 2 of the Galois group Gal(Ki/F) = (o"i) Il25l 

On the other hand, we have the cyclotomic algebra CA = (K 2 /L 2 ,cx 2 ■ Ci6 | — _ Ci6,72 = 
1 + Cs) = K 2 © u 2 K 2 , and cx 2 : Cie ^ Ki6,Cs ^ ~Cs- As we nave °~2(Cs) = -Cs, L 2 is the 
subfield fixed by (of) the subgroup of order [L 2 : F] = 2 of the Galois group Gal(K 2 /F) = (er 2 ) 
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From the nested sequence of fields (71), we can deduce that 

V 2 = -iV Kl/Ll (ai) -iJV Kl /Li(6i) - Cs Tr Kl/Ll (ai6i) . 




(65) 



e Li c Ki 



e Fi c K 

On the other hand, we have V\ G K, we can define it as Pi = a 2 + + Ci6 c 2 + Ci6^2 with 
a 2 , b 2 , C2, c?2 G Ki, then 



Pi + <t 2 (Pi)= 2a 2 +Ci6(l + i)62 + Ci 6 (l-*X 



(66) 



Therefore, we can define A^4 X 4 as 
^4x4 = V 1 + a 2 (V 1 )+V 2 

= {-N^/uiai) ~ iN Kl/u (bi) + 2a 2 ) + Cie(l + i)h - Cie^Ki/Li (oi&i) + Ciel 1 _ 0^2 
= A + Ci 6 5 + Ci 6 C + &Z> (67) 

with 



--Wki/Li(gi^ 



^1^(60 + 202 G Ki 



A = 

B = (l+i)b 2 G Ki 

C = Tr Kl / Ll (ai&i) G Li C Ki 

D = (1 - i)d 2 g Ki (68) 

It can be seen as a vector space of Ki with basis (1, Ci6, Ci6> Cie)» an( ^ tnus -^4x4 = if and 
only ifA = £>=C = .D = 0. This condition reduces to 

A = N Kl/hl (ax) + iN Kl/hl (bi) = 2a 2 (i) 

B = 6 2 = (ii) 

= (69) 
C = =4- Tr Kl/Ll (ai&i) = (in) 

D = rf 2 = (iv) 

So in order to prove that the 4x4 minor is non-zero, we have to prove that the latter condition 

cannot be verified. We proceed by contradiction. 

For this task, we show that by assuming that (ii), (Hi), (iv) are verified, we cannot have (i). 

In fact, if A = 0, one particular case would be when a\ = b\ = 0, so that a 2 = 0. 

However, if ax = b\ = and according to Equations (1631) and (in), we have 7-2 = Tr Kl / Fl (F) = 

0. Consider the general case where X G K, we can define it by 

X = a + Cie/3 with a, /3 G F 1? (70) 
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and its conjugate by a\{X) = a + CwP- We have then 

Y = XalaxiX) = (a + da/3) (a - CieyS) = (<*<* ~ ClM + Cie(a/3 - Cie/^)- (71) 
Since Tr K / Fl (F) = 0, thus Y G K is of the form Y = 7(1$, with 7 G F x . Therefore, we have 

act - Ci 6 /3/3 = =^ aa = Cs/3/3- (72) 

Let us now compute Pi + 02 (Pi) given this condition according to P2 = 0. Recall that 
Pi = -iM\o\o\{M\) with Mi = iXa 2 cr 1 (X)-a 1 (X)a 2 (X), then Mi and Pi can be reduced 
to 

Mi = -2aa - 2Ci 6 a/3 = -2a(a + CwP) (73) 

Pi = -4aa(a + Ci6^)(a-iCi 6 |) (74) 

and 

Pi+0 2 (Pi) = -Ma3>(2aa + - i)a/3 + (i 6 (l + 1)8*^. (75) 

On the other hand, we have according to (ii), (iv) (62 = ^2 = 0) and Equation (1661) that 
Pi + 02 (Pi) = 2a2- So, it can be simplified to 

Pi + 02 (Pi) = Siaaad. (76) 

Therefore, 

Pi + 0-2 (Pi) = a = 0. (77) 

However, a = means that .Mi = as well. But, we have already proved in Proof 2 that 
Mi 7^ over Z[i] unless s = 0. Consequently, Pi + 02 (Pi) 7^ 0, then a 2 7^ 0. So Given 
Oi = bi — 0, we prove here that A 7^ and thus, ^4x4 cannot be zero over Z[i] for s ^0. 

This last proof concludes the analysis on the full-rank asynchronous codeword matrices for 
the different types of delay profiles i.e., the full-diversity of the 4x4 code T, hence its delay- 
tolerance for arbitrary delay profiles. 

VI. New Delay-Tolerant Codes from other Perfect Codes 

We derive now delay-tolerant codes from the perfect codes presented in ifTTl . These latter 
codes differ from the previous ones by the construction of their generator matrices and their 
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non-norm element 71. Whereas this element was chosen as a root of unity in 2, 3, 4-dimensional 
perfect codes (7! = for the current codes it is of the form 

7T 

7i = — 

71* 

where ir is an element of Ki and n* its complex conjugate, ir is chosen as a suitable prime in 
Z[i] or Z[j] so that the element 7x is of unit norm and it is non-norm for the extension K^/F. 

Based on the same approach in Section (IIV-AI) . the delay-tolerant code is constructed using 
the tensor product of two number fields with the same degree and coprime discriminants. In 
previous dimensions, the second field corresponds to the cyclotomic extension K 2 = Q(C M ) 
where ( M is the M th root of unity since the non-norm element of the perfect code is itself a 
root of unity. Consequently, the relative extension will be here K 2 = Q(0 2 ) with 9 2 = Wji is 
the M th root of the non-norm element 71. 



A. 2 x 2 Code 



We consider the case of 2 antennas. The corresponding 2x2 perfect code was constructed 

in [HI on the field F = Q(i), and thus transmits g-QAM symbols. Let 6 X = 2cos(y L ) and 

Ki = Q(i, 6 X ) the relative extension of F of degree [Q(i, 6x) : Q(i)] = 2. The cyclic group Ki/F 

is generated by a\ and the cyclic algebra is then At(Ki/F, a±, 71) with the non-norm element 

3 + 2i 
71 ~ 2 + 3i 

The rotated lattice Z[i] 2 is obtained by a technique presented in IfTTll . E4ll different from the one 
used for previous perfect codes. The generator matrix is numerically given by Il26l 

-0.52573 -0.85065 
-0.85065 0.52573 

Now, let K 2 = Q(9 2 ) the cyclotomic extension of degree 2 of F with 62 = (71) 1 / 2 . Its relative 
discriminant is g?k 2 = 52 and is coprime to = 5. The cyclic Galois group generator is 
cr 2 : 6> 2 —8 2 and the integral basis is B 2 = (1,6*2)- The rotated Z[i] 2 lattice is generated by 
the unitary matrix 



Mi 



(78) 



M2 = 7I 



02 
-9 2 



(79) 



Then, the compositum of both cyclic extensions is defined by K = KiK 2 = Q(z, 6>i, (71 ) 1 / 5 ) 
of order 4 over Q(i) and accordingly the 4-dimensional complex lattice is generated by the 
4x4 unitary matrix M 4 = M 2 <g) M x . The codeword components are derived from the linear 
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combination x = M 4 s of g-QAM information symbols. They are then distributed in the 2x2 
codeword matrix assigned by the Fourier matrix coefficients (pi, I = 1, . . . , 4 in dimension 2 in 
order to guarantee the same NVD as the corresponding 2x2 perfect code. Both matrices U 
and V can also be derived by replacing ( 8 by 6> 2 . Moreover, the code construction allows to 
have a non-zero product distance, yielding a delay-tolerant code that maintains its full diversity 
regardless of the timing offset among its rows as shown for the previous 2x2 delay-tolerant 
code. 

VII. Performance Evaluation of Delay- Tolerant codes 

In this section, we evaluate the performance of the proposed distributed space-time codes used 
by the relays in synchronous as well as asynchronous transmission. Recalling the cooperative 
system model presented in Section IIII-B L a virtual MIMO scheme is assumed with M transmit 
antennas (one per relay) and N r receive antennas. The decoding is performed using the Sphere 
Decoder as for the perfect codes in conventional MIMO transmission. However in the case of 
asynchronous relays, the codewords are transmitted over T + D max symbol intervals resulting 
in rank deficiency of the channel matrix. In order to tackle this problem, the MMSE-DFE 
preprocessing ll27Tl is required to precede the lattice decoding so that the transformed channel 
has always full rank. 

The performance are represented in terms of codeword error rate CER and bit error rate BER 
versus signal-to-noise ratio Et/N per receive antenna, which is adjusted as 



dB 



iVn 



dB 



10 log R 



(80) 



where E s is the average energy per receive antenna and R is the code rate in bits per channel 
use (bpcu). 

A. Performance Comparison of Existent 2x2 Codes 

For 2x2 schemes, we consider the full-rate full-diversity existent space-time codes in this 
dimension, namely the Golden code G [14], or its variation matrix C proposed in [28] 



C(s) 



1 



Si + irs 4 rs 2 + s 3 
s 2 - rs 3 irs\ + s 4 



/2(l + r 2 ; 

the Silver code (Tirkkonen-Hottinen Code) ||29l , 11301 defined by 

X = X A ( ai ,s 2 ) + TWX fl (s 3 ,s 4 ) 



01-1, 



with ~K(si,Sj) 



Si Sj 




1 


, T = 








-1 



w 



1 

V7 



1+i -l + 2i 
1 + 2i l-i 
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the Sezginer-Sari code OTTl defined by 

S(s) 

the Damen code D [12] defined by 

D(s) 
with a 



as\ + bs3 —cs 2 — ds\ 
as 2 + bsi cs\ + dsl 

(l-y/7)+i(l + y/7) 

4%/2 



d = —ib. 



X 1 


sa 




asi + bs 2 - 


cs 3 - 


ds± 


—csi — ds 2 


— as 3 — bs^ 


x 2 


X4 




—bs\ + as 2 - 


h ds 3 


— CS4 


—dsi + cs 2 


— bsz + as4 



d 



and the 



a/(5+v / 5)(2+v / 2)' x /(5- v / 5)(2+V2)' • n /(5+n/5)(2- n /2) ' a/(5-v/5)(2-n/2) 

new proposed code T given in equation (fT9l ). These codes are compared in a distributed setup 
with and without delays. Note that the code D has been proved to verify the NVD criterion for 
any constellation carved from Z[i] and to be delay-tolerant It32l . 

In the above 2x2 schemes, the codewords matrices contain 4 modulated information symbols 
carved from 4-QAM constellation and transmitted over T = 2 channel uses. The transmission 
rate is hence R = ttt! — , where max = 0, 1 is the maximum delay with = (1,0) the delay 
profile in asynchronous transmission. 

Figure |4] shows the codes performances for synchronous relays (o max = 0). Observe that the 
Golden code G (or C) outperforms all the other codes. For example, it has about 1 dB and 0.5 
dB gains over D and T, S at a BER of 10 -4 , respectively. Note also that the new code T gives 
the same performance of the Golden code. 

Whereas for asynchronous relays, the situation is reversed between codes D and G for a delay 
of one symbol period since the latter is not delay-tolerant. It can be seen in Figure \5\ that both 
delay-tolerant codes T and D provide gains of 2 dB and more than 3 dB over codes T and G, 
S at a BER of 10~ 4 , respectively. In addition, it can be noticed that T performs almost similar 
to D, and it merely improves for high SNR (> 13 dB) by 0.2 dB at a BER of 2 x 10~ 5 . 

Using the unitary matrices U and V that provide the new code T from the Golden code G, 
we can also obtain new delay-tolerant codes based on T and S codes as 



UTV and S d = USV 



(81) 



Note that U and V are not necessarily the optimal matrices for these codes, but they allow to 
have new delay-tolerant codes with the same determinants as the initial ones. One can easily 
verify as demonstrated for code T that the product distances associated with these new codes 
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are non-zero over Z[i]. Figure [6] depicts the performances of the new codes for asynchronous 
relays with a delay of 1 symbol period. It can be observed that all these delay-tolerant codes 
preserve their diversity and that the code T gives the best performance. At a BER of 2 x 10~ 5 , 
it gains about 0.2 dB and 0.8 dB over and S^, respectively. 

B. Performance o/3 x 3 Codes 

For the 3x3 schemes, 9 modulated symbols carved from 4-HEX constellation (e 
are transmitted at a rate of R = ^— ^ — bpcu, where D max = 0, 2 is the maximum delay and 
D = (2, 1, 0) the delay profile in asynchronous transmission. 

In Figure |7J we can observe that both the perfect code and the new code T have the same 
performance for synchronous relays. Whereas for asynchronous relays, the delay-tolerant code 
preserves the diversity and provides a gain of 5 dB over the 3x3 perfect code at CER of 10~ 4 
for max = 2. 

VIII. Conclusion 

In this paper, we have proposed new delay-tolerant space-time codes based on the perfect codes 
algebras. Using tensor product of the perfect code field extension with another field extension 
of the same order M over the same base field and which Galois extensions have coprime 
discriminants, we build rotated lattices in higher dimension in order to construct MxM codes. A 
key parameter in the construction is the coefficients (pi that allow to preserve the same properties 
of the perfect codes in synchronous transmission. 

We have found that <pi corresponding to the coefficients of the Fourier matrix in dimension 
M yield the same non-vanishing determinants as the perfect codes. These codes besides having 
full-rate, full-diversity, uniform energy per transmit antennas {\<pi\ = 1) and are information 
lossless, they have the NVD property and thus are optimal DMT achieving in synchronous case. 

In addition, for asynchronous transmission, we have proved for M = 2, 3, 4 that the new codes 
preserve their full-diversity and are delay-tolerant for arbitrary delay profiles. This property is 
obtained thanks to the non-zero product distances over Z[i] or Z[j] and the full-rank minors of 
the delayed matrices. 
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Virtual M x N r MIMO 



Fig. 1. Cooperative System Architecture 




Fig. 2. Compositum Field of the Tensor Product Algebra 
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Fig. 3. Nested Sequences of fields 




Fig. 4. Performances of the codes T, D, G (or C), T and S without delay 
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Fig. 5. Performances of the codes T, D, G (or C), T and S with a delay of 1 symbol period 




Fig. 6. Performances of the codes T, D, Td and S<j with a delay of 1 symbol period 
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Fig. 7. Performances of 3 x 3 codes w/wo delay 



